Cell Metabolism

Intermuscular adipose tissue (IMAT) expansion is closely associated with cardiometabolic disease, yet its cellular organization and regulatory mechanisms remain poorly defined. Here, we define a human IMAT gene signature using bulk transcriptomics and identify candidate regulators for IMAT function including adipogenic transcription factor early B-cell factor 2 (EBF2). To determine how these programs are organized in situ, we mapped this signature in a mouse model of diet-induced CMD using spatial transcriptomics. We found that IMAT expansion occurs within discrete stromal niches surrounding muscle fibers, characterized by coordinated activation of adipogenic, extracellular matrix, inflammatory, and metabolic pathways. Spatial analyses showed that fibro-adipogenic progenitor (FAP) abundance does not predict adipocyte formation, supporting a model of localized and context-dependent lineage transitions. Cross-species comparison revealed partial conservation of human IMAT gene programs, validating the mouse model and highlighting species-specific features. Functional experiments in human primary myoblasts showed that EBF2 is sufficient to induce adipogenic reprogramming. Our findings establish IMAT as an active, spatially organized remodeling niche and identify lineage plasticity as a central mechanism driving its expansion in metabolic disease

GLP-1 receptor agonists induce substantial weight loss, but the extent to which lean tissue and physical function are preserved in routine care remains poorly understood. Using an EHR-linked body-composition digital phenotyping pipeline with LLM-based extraction, we performed a large-scale longitudinal analysis of 670,422 first-episode GLP-1RA users, including 456,742 treated with semaglutide and 213,680 treated with tirzepatide. Among these, 7,965 individuals with paired pre- and post-initiation body-composition measurements were analyzed over 12 months. Tirzepatide was associated with greater relative lean body mass (LBM) loss than semaglutide at each measured time point, with excess LBM losses of 1.1%, 1.5%, 1.3% and 2% at 3, 6, 9 and 12 months, respectively. A Depletive GLP-1 metabotype, defined as >20% total body weight (TBW) loss with >5% LBM loss, was significantly more frequent with tirzepatide than semaglutide during the first year of therapy (10.3% versus 6.7%, p<0.001). By contrast, a Prime GLP-1 metabotype, defined as >10% TBW loss with <5% LBM loss, was numerically more frequent with semaglutide than tirzepatide, but not significantly so (12.3% versus 11.8%, p=0.66). Higher drug dose and longer exposure were associated with progressively greater LBM decline in both treatment groups (both p<0.001). Among 3,746 examined EHR phenotypes, baseline musculoskeletal pain emerged as the most significant correlate of greater LBM loss (BH-adjusted q<0.001): cervicalgia (semaglutide, -4.1 percentage points; tirzepatide, -14.3 percentage points) and knee pain (semaglutide, -4.8 percentage points; tirzepatide, -13.4 percentage points), consistent with mobility-limited patients being more vulnerable to lean-tissue depletion during incretin therapy. Analysis of EHR notes for on-treatment functional features showed reduced exercise tolerance was the strongest correlate of greater LBM loss, increasing by 7.2 and 11.1 percentage points in semaglutide- and tirzepatide-treated patients, respectively. An independent analysis of all available Single-cell RNA-seq data from human musculature showed broader GIPR+ cellular distribution than GLP1R+ cells across immune, stromal, vascular, and contractile compartments, providing plausible biological context for the greater LBM loss observed in routine care with tirzepatide (dual GLP1R-GIPR agonist) relative to semaglutide (GLP1R-specific agonist). In this observational study, greater weight-loss efficacy did not necessarily translate into more favorable body-composition outcomes, underscoring the need for clinical decision-making and trial designs that maximize each patient's likelihood of achieving a Prime GLP-1 metabotype.

The transcription coregulator OCA-B promotes CD4+ T cell memory recall responses and autoimmunity. OCA-B T cell deletion prevents spontaneous type-1 diabetes (T1D) onset in non-obese diabetic (NOD) mice and blunts T1D in a subset of more aggressive models. However, the role of OCA-B in diabetes induced by treatment with immune checkpoint inhibitors (ICIs), and the role of OCA-B in the control of tumors with and without ICI treatment, has not been studied. Here we show that islet and pancreatic lymph node T cells from T1D individuals express measurable POU2AF1 mRNA. Deletion of OCA-B in T cells fully insulates 8-week-old non-obese diabetic (NOD) mice against ICI-induced diabetes and partially protects 12-week-old mice. Salivary and lacrimal gland infiltration and inflammation were also reduced. Protection was associated with a block in the differentiation of progenitor exhausted CD8+ T cells (TPEX) into terminally exhausted CD8+ T cells (TEX). We show that OCA-B T cell loss preserves anti-tumor immune responses following PD-1 blockade in different tumors and mouse strains. These findings point to a potential therapeutic window in which pharmaceuticals targeting OCA-B could be used to block the emergence of both spontaneous and ICI-induced autoimmunity while sparing anti-tumor immunity. We develop first-in-class small molecule inhibitors of Oct1/OCA-B transcription complexes and show that administration into NOD mice also blocks diabetes emergence following PD-1 blockade. These results identify OCA-B as a promising therapeutic target for the prevention of autoimmunity and immune-related adverse events (irAEs).

Dietary patterns worldwide have shifted toward increased consumption of ultraprocessed foods (UPFs), which has been linked to higher disease burden. One mechanism proposed to impact both their consumption and contribution to metabolic disease is altered post-ingestive metabolic response in comparison to nutritionally similar foods. Here, we recruited 57 healthy-weight 18-45-year-old adults to examine the effects of food processing on postprandial metabolism and brain response. Despite nutritional matching, UPF meals evoked a greater insulinemic and energetic response with attenuated carbohydrate oxidation relative to non-UPF meals. Next, between-condition differences in peak carbohydrate oxidation were associated with mesolimbic and superior temporal gyrus activation in response to food cues. Finally, although food value did not differ between conditions, brain responses correlated with food valuation were positive for non-UPF but negative for UPF in visual cortex and striatum. These findings demonstrate that food processing influences post-ingestive metabolism in a way that could help explain long term health effects and differences in food reward through mechanisms beyond calories and macronutrient composition alone.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the deadliest human cancers. The current largest published PDAC Genome-Wide Association Study (GWAS) identified 23 genetic risk signals, but most lack sufficient characterization. This study aimed to functionally characterize the chr13q12.2 (PLUT/PDX1) PDAC GWAS risk locus. Fine-mapping, luciferase reporter assays, and electrophoretic mobility shift assays implicated rs9581943, a PDX1 promoter SNP, as a functional variant underlying this GWAS signal. GTEx expression QTL analyses identified rs9581943 as a significant PDX1 eQTL in pancreas, and CRISPR/Cas9 editing in PDAC-derived cell lines confirmed a functional relationship. PDX1 is a transcription factor involved in early pancreas development and {beta}-cell homeostasis, but its role in exocrine pancreatic cells is unclear. Single-nucleus RNA-seq analyses of pancreatic acinar and ductal cells from neonatal, adult, and chronic pancreatitis donors suggested PDX1 activity alleviates high secretory load and ER-stress in acinar and biases ducts toward homeostatic phenotypes. Similarly, scRNA-seq analyses of pancreatic tumors suggested PDX1 activity reduces biosynthetic and inflammatory stress and promotes epithelial differentiation. Our study therefore implicates rs9581943 as a causal variant for the chr13q12.2 PDAC GWAS signal wherein the risk allele reduces PDX1 expression, eroding PDX1's capacity to buffer stress and stabilize epithelial cell fate in the exocrine compartment.

Semaglutide has shown benefit in metabolic dysfunction-associated steatohepatitis (MASH), but real-world evidence across longitudinal liver phenotypes remains limited, particularly regarding how liver remodeling relates to weight loss and dose exposure. Using a de-identified federated electronic health record network spanning more than 29 million patients in the United States, including 489,785 semaglutide-treated adults, we analyzed 6,734 patients with baseline liver disease burden. We find that higher attained pre-landmark (0-2 years) semaglutide dose was associated with lower post-landmark (2-4 years) risk of steatohepatitis, alcoholic liver disease, and all-cause mortality, whereas greater pre-landmark weight loss was associated with lower post-landmark risk of steatohepatitis, steatotic liver disease, and hepatorenal syndrome, indicating distinct dose- and weight-linked patterns of long-term liver benefits. These associations were notable because semaglutide prescribing was generally lower during the post-landmark period, raising the possibility of durable benefit beyond peak exposure. Towards better understanding mechanistic bases for liver protection, we performed a complementary longitudinal study of 326 adults with paired noninvasive liver elastography measurements before and after treatment initiation. Median liver stiffness decreased from 4.85 [3.02 - 7.20] to 3.9 [2.6 - 5.8] kPa after semaglutide initiation (median change = -0.38 kPa; p<0.001), with 194 of 326 patients (59.5%) showing lower follow-up stiffness. A clinically meaningful reduction of at least 20% was observed in 133 of 326 patients (40.8%), and 69 of 326 (21.2%) shifted to a lower fibrosis stage by prespecified elastography thresholds. Larger improvements were also seen in patients with higher baseline stiffness (p<0.001); notably 80% of patients with cirrhosis-range baseline stiffness ([&ge;]12.5 kPa) achieved [&ge;]20% improvement versus 29.5% with minimal baseline disease (p <0.001). The proportion achieving at least 20% stiffness improvement was similar across weight-loss strata, including patients with no weight loss or weight gain and those with at least 10% weight loss (38.0% in each group), and liver stiffness change showed negligible correlation with changes in weight, BMI, HBA1c, alanine aminotransferase, or aspartate aminotransferase. To provide biological context, single cell RNA analyses demonstrated sparse overall hepatic GLP1R expression (0.0239%), with enrichment in non-parenchymal niches including cholangiocytes, intrahepatic cholangiocytes, liver sinusoidal endothelial cells, and hepatic stellate cells implicated in fibrogenesis and vascular remodeling. Together, this real-world evidence suggests diverse liver benefits for semaglutide beyond weight-loss with intricate dose response relationships.

Dementia incidence has been declining in Western societies for decades, but whether this reflects higher cognitive capacity entering old age, slower cognitive decline, or both remains unresolved. Analysing ~783,000 episodic memory assessments from ~219,000 individuals across five longitudinal cohorts, we find that later-born cohorts benefit from a double dividend: higher memory levels entering old age and slower rates of decline. The projected 20-year cohort advantage at age 80 is of sufficient magnitude to plausibly account for the observed 13% per-decade decline in dementia incidence reported in meta-analyses. Generational gains are disproportionately concentrated among the fastest-declining individuals, and are reflected in lower hippocampal atrophy rates in an independent sample. A formal bounding analysis shows that the double dividend is robust across a range of plausible period assumptions, consistent with environmental conditions operating across the lifespan having reshaped the architecture of human cognitive aging.

Amyotrophic lateral sclerosis (ALS) is a rapidly progressing neurodegenerative disease with a heterogeneous clinical presentation, complicating early diagnosis and therapeutic monitoring. To identify disease-specific biomarkers, we performed an unbiased cerebrospinal fluid (CSF) proteomic analysis in 87 ALS patients, 89 healthy controls, and 61 neurological controls using mass spectrometry. Across all quantified proteins, 399 were significantly dysregulated in ALS, including established neurodegeneration (NEFL, NEFM, UCHL1) and neuroinflammatory (CHIT1, CHI3L1, CHI3L2) markers. Correlation and pathway analyses uncovered dysregulation of immune, synaptic, and metabolic processes, with aberrant complement activation emerging as a hallmark. Complement proteins increased progressively with declining ALS Functional Rating Scale-Revised and longer disease duration, whereas early-stage markers (CLSTN3, CHAD, RELN) indicated pre symptomatic neuronal and synaptic disruptions. Machine learning identified a minimal five protein CSF panel (MB, ITLN1, YWHAG, FCGR3A, PGAM1) that accurately distinguished ALS patients from healthy controls, capturing disease-specific pathophysiology beyond general neurodegeneration. Our findings define a robust ALS-specific CSF proteomic signature, reveal prognostic protein candidates across disease stages, and provide a framework for diagnostic biomarker development, enabling earlier intervention and monitoring.

Quantitative measures for tracking functional health have generally been lacking. Intrinsic capacity (IC) has been proposed as an appropriate measure, but its metrics have been derived in small datasets and sparse longitudinal data. Using harmonized measures of cognition, locomotion, sensory function, vitality, and psychological well-being from 501,615 UK Biobank participants and followed for a median of 15.5 years, we derived domain-specific and composite IC scores. We examined associations with incident disease, cause-specific mortality, multimorbidity, lifestyle and socioeconomic factors, and multi-omic profiles from Olink proteomics, NMR metabolomics, clinical biochemistry, and blood-cell traits. We found that composite IC declined non-linearly with age, and within-person decline was steeper than the cross-sectional age measures. Participants with greater baseline morbidity, those who subsequently developed incident disease, and those who died earlier in follow-up showed lower IC trajectories across adulthood. The IC domains were only modestly correlated with one another, supporting multidimensionality, yet higher overall IC was associated with lower risk of most diseases examined. The dominant IC domain varied by endpoint, with cognition informative for dementia, sensory function for hearing loss, psychological capacity for depression, locomotion for osteoarthritis, and vitality for cardiometabolic outcomes. IC was also associated cross-sectionally with physical activity, insomnia, smoking, medication burden, and socioeconomic disadvantage. More proteins were found predictive for vitality, and enrichment converged on immune/inflammatory and metabolic pathways. Blood-based surrogates recapitulated part of the phenotypic signal, particularly for vitality. Overall, this IC framework captures longitudinal health trajectories and broad disease vulnerability in a large middle- to older-aged cohort and supports IC as a clinically meaningful, multidomain phenotype of aging and identifies blood-based correlates that may facilitate at-scale future monitoring of aging-related function declines.

Background and Aims: Alterations in immunoglobulin G (IgG) N-glycosylation are implicated in inflammatory bowel disease (IBD); however, the robustness of IgG glycan signatures across IBD cohorts with diverse demographics and geographic origins remains underexplored. We aimed to determine whether compositional data analysis (CoDA) and machine learning (ML) can identify IBD-related IgG N-glycan signatures and whether these signatures capture disease-associated acceleration of biological aging. Methods: We analyzed the IgG glycome profiles of 1,367 plasma samples collected from healthy controls (HC), symptomatic controls (SC), and people with newly diagnosed Crohn's (CD), and ulcerative colitis (UC) across four cohorts (UK, Italy, United States, and Netherlands). IgG glycosylation was analyzed by ultra-high-performance liquid chromatography, yielding 24 total-area-normalized glycan peaks (GPs). Analyses were performed using cross-sectional data obtained at baseline. CoDA-powered association analyses were used to identify disease-related effects on GPs while controlling for demographic covariates. ML models were trained and evaluated to assess generalizability to unseen cohorts and demographic subgroups, with a focus on discrimination and reliability. Results: Across all cohorts, people with IBD demonstrated accelerated biological aging as quantified by the GlycanAge index. This was accompanied by consistent reductions in IgG galactosylation, with effects partially modulated by age. Classification models trained on glycomics and demographics achieved robust discrimination (AUROC~0.80) between non-IBD (HC+SC) and IBD across cohorts. Conclusion: These findings reveal accelerated biological aging in people with IBD and support the translational potential of IgG glycans as biomarkers and a novel route toward clinically interpretable personalized risk estimates.

Background: The glucagon-like peptide-1 receptor (GLP1R) is a key regulator of glucose metabolism and appetite and a major therapeutic target for type 2 diabetes (T2D) and obesity. Genetic studies have implicated the GLP1R locus in both body mass index (BMI) and T2D, but it remains unclear whether their underlying genetic associations are the same. Methods: We analyzed 431,107 participants of genetically inferred European ancestry from the Million Veteran Program. Within 500 kb of GLP1R, we performed locus-wide linear regression models for BMI and logistic regression models for T2D, adjusted for age, sex, and 10 principal components. We identified primary and secondary BMI sentinel variants using conditional analyses and evaluated their associations with T2D. Bayesian fine-mapping was used to construct credible sets of GLP1R locus for BMI and T2D. Results: Conditioning on the primary sentinel variant rs12213929 (upstream of GLP1R, {beta} = 0.11; 95% CI 0.09-0.14; p = 1.94E-17), we identified a secondary variant (rs13216992, intron of GLP1R) independently associated with BMI ({beta} = 0.10; 95% CI 0.07-0.13; p = 7.88E-14). The two sentinel variants showed low linkage disequilibrium (r2 = 0.03). A two-variant allelic burden score (0-4; sum of the rs12213929 G-allele count and rs13216992 C-allele count) showed that participants with 4 risk alleles had 0.47 kg/m2 higher BMI than those with 0 risk alleles (95% CI 0.39-0.55; p < 2E-16). Both variants were associated with higher T2D risk, but with distinct patterns after BMI adjustment: the rs12213929-T2D association persisted after adjustment for BMI (OR = 1.02; 95% CI 1.01-1.03; p = 0.0004), whereas the rs13216992-T2D association was fully attenuated (OR = 1.00; 95% CI 0.99-1.01; p = 0.68). Fine-mapping identified a compact 95% BMI credible set of 17 variants and a broader 95% T2D credible set of 42 variants, with all BMI credible variants contained within the T2D set. Conclusions: The GLP1R locus harbors at least two independent BMI-associated variants that exhibit heterogeneous relationships with T2D: rs12213929 influences T2D risk partly through BMI-independent pathways, whereas rs13216992 appears to act predominantly via adiposity. These findings refine the genetic architecture at this key therapeutic target gene and provide a foundation for functional and pharmacogenomic studies to determine whether GLP1R variation can inform precision prevention and treatment of obesity and T2D.

Background. Detection of cerebral amyloid pathology currently requires amyloid PET imaging ($5,000-$8,000) or cerebrospinal fluid analysis via lumbar puncture, procedures that are inaccessible for population-level screening. The FDA-cleared Lumipulse G pTau217/Abeta1-42 plasma ratio test (May 2025) represents the first approved blood-based alternative; however, single-ratio approaches cannot distinguish Alzheimer's disease (AD) from non-AD neurodegeneration or provide multi-dimensional disease characterization. Methods. We developed Virtual Spectral Decomposition (VSD), a framework that decomposes plasma biomarker profiles into biologically interpretable diagnostic channels. Four plasma biomarkers - phosphorylated tau-217 (pTau217), amyloid-beta42/40 ratio, neurofilament light chain (NfL), and glial fibrillary acidic protein (GFAP) - were measured in 1,139 Alzheimer's Disease Neuroimaging Initiative (ADNI) participants. Each biomarker was mapped to a VSD channel representing a distinct pathophysiological axis: tau/amyloid phosphorylation, amyloid clearance, neurodegeneration, and astrocytic activation. Channel weights were calibrated via logistic regression, and performance was evaluated against amyloid PET (UC Berkeley) using 10x5-fold repeated cross-validation. Results. VSD 4-channel fusion achieved AUC = 0.900 (+/-0.018), exceeding pTau217 alone (0.888+/-0.022). Optimal sensitivity was 89.7% with 78.1% specificity (NPV = 90.8%). The NfL channel received a negative weight (beta = -1.1), functioning as a disease-exclusion signal: elevated neurodegeneration without amyloid-tau coupling actively reduces the AD probability, distinguishing AD from non-AD neurodegeneration. Complementary CSF proteomics analysis (7,008 proteins, 533 participants) identified 17 amyloid-specific proteins (0.24% of the proteome), revealing a 49:1 tau-to-amyloid asymmetry that explains why blood-based tau markers outperform amyloid markers. Conclusions. Blood-based VSD provides an interpretable, multi-channel framework for amyloid detection that incorporates explicit disease-exclusion logic unavailable to single-biomarker approaches. The architecture extends to multi-disease screening, where the same blood specimen could be routed through disease-specific modules for AD, Parkinson's disease, and cancer.

Colorectal cancer (CRC) is the second leading cause of cancer death globally and the number one cause of cancer death in people under 50 years old. The reasons for the rise of early-onset CRC are unknown, and while anatomically distinct subtypes of CRC have substantial clinical and molecular associations, the etiology of region-specific disease, such as early-onset CRC's enrichment in the distal colon, remains unclear. Understanding regional mutagenesis may identify risk factors for this public health concern and CRC more broadly. To evaluate mutational dynamics across the premalignant colon, we performed whole-genome sequencing of 125 individual colon crypts taken from six standardized regions biopsied during colonoscopy, collected from 11 donors without polyps and 10 with polyps. We observed mutation spectra and accumulation rates consistent with previous whole-organ studies, with greater subclonal mutation capture enabled by experimental design. T>[A,C,G] mutations, which are associated with colibactin genotoxicity from pks+ Escherichia coli, were significantly enriched in the rectum of donors with and without polyps (adjusted p-values < 0.01). Moreover, when comparing findings to crypts from individuals with CRC and sequenced CRC tumors, we observed consistent enrichment of the colibactin-associated mutational signature "ID18" in the rectum in both normal colon crypts and CRC tumors, without significant difference in colibactin-specific single nucleotide variant or insertion-deletion burden in crypts across the three clinical groups (i.e., no polyp, polyp, and CRC). These findings argue against a causal or prognostic role for colibactin in CRC, instead indicating that the proposed association with early-onset disease reflects anatomic specificity rather than cancer-specific clinical relevance.

Drug-tolerant persisters (DTPs) represent a major obstacle to durable responses in targeted cancer therapy. DTPs are commonly described as distinct single-cell states that survive drug treatment via reversible, non-genetic mechanisms and drive tumor recurrence. Recent work demonstrates that multiple DTPs can coexist, reflecting diversity in lineage, signaling programs, or stress responses. However, each DTP is still generally viewed as a uniform cellular phenotype. Building on our prior work describing a population-level DTP termed "idling" [Paudel et al., Biophys. J. (2018) 114, 1499-1511], here we present evidence supporting a fundamentally different view: that DTPs are not single-cell states, but rather heterogeneous populations composed of multiple sub-states with distinct division and death rates that balance to produce near-zero net population growth. Using single-cell transcriptomics and lineage barcoding, we identify multiple phenotypic states within idling DTP populations, with reduced heterogeneity compared to untreated populations, and find that idling DTP cells emerge from nearly all lineages. Transcriptomic and functional analyses further reveal altered ion-channel activity in idling DTPs, which we confirm experimentally. Moreover, drug-response assays reveal increased susceptibility of idling DTPs to ferroptosis, a non-apoptotic form of regulated cell death, indicating the emergence of vulnerabilities associated with drug tolerance. Altogether, our results support a population-level view of tumor drug tolerance in which DTPs comprise stable collections of phenotypic states, shaped by treatment-defined phenotypic landscapes, which are potentially vulnerable to subsequent interventions. This perspective implies that eradicating DTPs will require a fundamental shift away from cell-type-centric strategies toward sequential treatments that progressively reduce phenotypic heterogeneity by modulating the molecular and cellular processes that establish the DTP landscape, an approach previously termed "targeted landscaping."

Pathogenic variants in TP53, the key tumour-suppressor gene underlying Li-Fraumeni syndrome (LFS), are among the best-established causes of inherited cancer predisposition. However, large-scale sequencing has revealed that many apparently pathogenic TP53 variants detected in blood are the result of somatic clonal expansions, complicating risk interpretation. Using blood-derived whole-exome data from 469,391 UK Biobank participants, we combined variant allele fraction (VAF) with haplotype-sharing analysis to distinguish germline and somatic TP53 variants. Germline variants were concentrated at sites linked to partial loss of p53 function and lower disease penetrance, whereas classic LFS alleles appeared almost entirely somatic. High-VAF carriers of classic LFS alleles conferred markedly increased risk of haematological malignancy but not solid tumours, consistent with large TP53-mutant clonal expansions. The prevalence of somatic clonal expansion also correlated with missense variant pathogenicity, suggesting that somatic activity provides an informative in vivo proxy for functional impact. These results provide new insights into TP53-associated cancer risk at the population level, demonstrate that somatic rather than germline risk predominates in middle-aged healthy adults and provide a scalable framework for variant classification in large-scale population genomics.

Effective predictive biomarkers for immune checkpoint inhibitor (ICI) therapy remain an unmet need across solid tumors. Here, we present an integrated spatial proteomics workflow that combines in situ proximity ligation assay with multiplexed immunofluorescence to directly resolve PD1/PDL1 signaling events at the level of defined cellular phenotypes and their spatial organization within intact tumor tissue. Applied as a proof of concept to tumor samples from patients with metastatic urothelial carcinoma treated with pembrolizumab, this approach reveals that PD1/PDL1 interactions specifically involving cytotoxic CD8CD3 T cells are significantly enriched in complete responders, while such interactions are rare in patients with progressive disease. This interaction defined T cell subset achieves superior discrimination of clinical response compared to single marker PDL1 expression or immune cell abundance alone. By integrating direct detection of protein protein interactions with high dimensional single cell phenotyping, our workflow provides a mechanistically informed, spatially resolved biomarker of functional immune engagement. Beyond urothelial carcinoma, this platform establishes a generalizable framework for translating spatial signaling biology into predictive tools for immunotherapy response across tumor types.

De novo mutations (DNMs) arising in the parental germline are a major cause of severe developmental disorders. While most DNMs originate in the paternal germline, it remains unclear whether fathers of affected children carry a systematically altered burden of transmissible germline risk, or whether disease largely reflects stochastic outcomes of shared population-wide mutational processes. Here, we combined whole-genome sequencing of 168 parent-child trios with ultra-accurate duplex sequencing of paternal sperm to directly relate transmitted DNMs to the broader mutational and selective landscape of the male germline. In 127 fathers, sperm mutation burden and mutational spectra were indistinguishable from population reference cohorts. Positive selection metrics were likewise concordant, with a global dN/dS of 1.56 (95% CI 1.45-1.67) compared to 1.44 (95% CI 1.17-1.77) in controls and 28 of 32 significantly selected genes overlapping with prior findings. Six fathers harboured a pathogenic early mosaic variant detectable in sperm at allele fractions that ranged from 0.7% to 14.8%. Although these variants generated substantial individual-level risk outliers, they accounted for only [~]11% of the aggregated exome pathogenic burden across the cohort. The remaining burden was distributed across low-VAF mutations, including positively selected driver variants and other rare mutations accumulating with paternal age. Together, these results show that transmissible de novo disease risk is governed primarily by universal germline mutational and selective processes, while early developmental mosaicism produces uncommon but clinically meaningful deviations. This integrated view clarifies how mutation timing, age-associated accumulation and germline selection jointly shape inheritance risk.

Genome-wide association studies (GWAS) have transformed our understanding of human biology, but are constrained by the need for predefined phenotypes. We introduce Vector2Variant (V2V), a general-purpose framework that transforms any set of high-dimensional measurements (such as machine learning embeddings) into a genome-wide scan for associations, without requiring rigid specification of a phenotype. Rather than testing genetic variants against single traits, V2V finds the axis in multivariate space along which carriers and non-carriers maximally differ, and produces a continuous "projection phenotype" that can be interpreted by association with disease labels. The projection phenotypes correlate with orthogonal clinical biomarkers never seen during training, suggesting the learned axes capture biologically meaningful variation. We applied V2V to imaging, timeseries, and omics modalities in the UK Biobank and recovered established biology (like the role of CASP9 in renal failure) without the need for targeted measurements, alongside novel associations including a frameshift variant in LRRIQ1 (potentially protective for cardiovascular disease). V2V is computationally efficient at genome-wide scale, producing summary statistics and disease associations that facilitate target prioritization without the need for phenotype engineering.

Osteoarthritis (OA) is a prevalent musculoskeletal disorder and a leading cause of global disability. Although meniscal damage is a major risk factor of OA pathogenesis, genetic regulatory studies have remained largely confined to articular cartilage. Here, we establish the first comprehensive expression quantitative trait locus (eQTL) map integrating whole-genome sequencing and bulk transcriptomics from human meniscus (n=112) and cartilage (n=113). Supported by single-nucleus multiomics (cartilage: 56,549 nuclei; meniscus: 34,343 nuclei), we uncovered highly tissue-specific genetic risk architectures. Colocalization with OA GWAS identified 27 meniscus-specific, 28 shared, and 20 cartilage-specific causal genes. Chromatin-informed fine-mapping and deconvolution elucidated distinct pathogenic mechanisms; notably, meniscus-specific signals converged on VEGFA via rare promoter variants and an enhancer in fibrochondrocyte progenitors, alongside a shared eQTL for CLEC18A. Exploratory analysis suggested candidate compounds to reverse pathogenic gene expression. Our findings underscore the meniscus as a distinct genetic driver, molecularly reinforcing OA as an entire joint organ failure.

How embryonic cells generate large clones of cells in the adult represents a fundamental question in biology. Here, using melanocyte stem cells (McSCs) in the zebrafish as a model, we explore the function of the master melanocyte transcription factor (MITF) in safeguarding McSCs in embryonic development and their potential to pigment large clones in the adult. MITF is well known is for its role in the specification of melanoblasts from the neural crest (NC) and their differentiation into melanocytes, yet little is known about how this activity shapes the stem cell lineages. Here, we use live imaging coupled with single-cell transcriptomics and lineage tracing to show that MITF (mitfa in zebrafish) protects the melanocyte stem cell (McSC) fate in zebrafish. Utilizing a temperature sensitive mitfavc7 mutant, we show loss of Mitfa leads to a surprising premature and aberrant expansion of McSC progeny at the niche during embryogenesis, coupled with novel emergent transcriptional cell states. Linage tracing of McSCs from the embryonic to juvenile stages reveals Mitfa activity is subsequently required in regeneration by Schwann cell-like and melanocyte stem cell progenitors that serve as a reservoir for fast-responding pigment progenitors. Thus, the impact of Mitfa loss on the melanocyte lineage is cell-state and stage-specific. The emergent cell states upon mitfa loss may have important implications for our understanding the loss of MITF activity in human genetic disease and melanoma.